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Verfahren zur treibstoffoptimalen Ansteuerung einer Dusenanordnung 

Die vorliegende Erfindung betrifft ein Verfahren zur treibstoffoptimalen Ansteuerung 
einer Anordnung von Dusen an einem Raumfahrzeug. 

Ein solches Verfahren ist beispielsweise aus US 6,347,262 B1 fur den Fall eines 
spinstabilisierten Raumfahrzeuges bekannt. Eine Anordnung von Dusen an einem 
Raumfahrzeug, wie sie die vorliegende Erfindung betrachtet, dient insbesondere zur 
Bahn- und Lageregelung des Raumfahrzeuges. Eine soiche Bahn- und Lagerege- 
lung durch Dusen ist beispielsweise aus EP 0 750 239 A2 bekannt. 

Aus der EP 0 977 687 ist ein spezielles Verfahren zur treibstoffoptimalen Ansteue- 
rung einer Anordnung von Dusen an einem Raumfahrzeug bekannt, wobei zum 
Auffinden einer treibstoffoptimalen Losung fur die Ansteuerung die Losung eines 
konvexen linearen Optimierungsproblems erfolgt mit 

- einer Initialisierungsphase zum Auffinden einer ersten zulassigen Losung des 
linearen Optimierungsproblems und 

- einer anschliefienden Iterationsphase, in der ausgehend von der zulassigen 
Losung des linearen Optimierungsproblems eine iterative Optimierung eines 
Gutekriteriums erfolgt. 

Es wird bei diesem Verfahren ein dualer Simplexalgorithmus angewendet, der durch 
ein weitgehend ungerichtetes Suchverfahren eine moglichst optimale Losung fur das 
Problem finden soli, wobei es aber bei diesem Verfahren vorkommen kann, dass es 
fur einen aktuell vorliegenden Krafte-Momenten-Vektor und die vorliegende Dusen- 
anordnung keine Losung gibt. 

Aus N.Karmakar: A new polynominal time algorithm for linear programming, 
Combinatorica 4 (4), 1984, S. 373 - 395 ist ein grundsatzliches Verfahren zur 
Losung linearer Optimierungsprobleme in allgemeiner Form bekannt. 
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Aufgabe der vorliegenden Erfindung ist es, ein Verfahren zur treibstoffoptimalen 
Ansteuerung einer Anordnung von Dusen an einem Raumfahrzeug bereitzustellen, 
das eine zielgerichtete Suche einer auf jeden Fall zulassigen Losung des linearen 
Optimierungsproblems erlaubt. Diese Aufgabe wird gelost durch die Merkmale des 
Patentanspruchs 1 . 

Die vorliegende Erfindung betrifft ein Verfahren zur treibstoffoptimalen Ansteuerung 
einer Anordnung von Dusen an einem Raumfahrzeug, wobei zum Auffinden einer 
treibstoffoptimalen Losung fur die Ansteuerung die Losung eines linearen 
Optimierungsproblems, insbesondere eines konvexen linearen Optimierungsprob- 
lems, erfolgt mit 

- einer Initialisierungsphase zum Auffinden einer ersten zulassigen Losung des 
linearen Optimierungsproblems und 

- einer anschliefienden Iterationsphase, in der ausgehend von der zulassigen 
Losung des linearen Optimierungsproblems eine iterative Optimierung eines 
Gutekriteriums erfolgt. 

Gemafi der Erfindung ist nun vorgesehen, dass 

- bei jedem Iterationsschritt die Bildung eines skalierten Iterationsgradienten 
erfolgt und 

- der Iterationsgradient mit einem Begrenzungsfaktor fur eine maximale Itera- 
tionsschrittweite multipliziert wird, welcher unter Berucksichtigung mindestens 
einer Randwertbedingung fur eine zulassige Losung gebildet wird. 

Durch die Anwendung eines skalierten Iterationsgradienten im Gegensatz zu einem 
blofien Suchverfahren wie beim Stand der Technik kann ein gerichtetes Auffinden 
der optimalen Losung erfolgen. Bei der Bildung des skalierten Iterationsgradienten 
kann insbesondere wiederum mindestens eine Randwertbedingung fur eine 
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zulassige Losung mit einflieflen. Es kann mit der Anwendung eines skalierten 
Iterationsgradienten auch weitgehend ausgeschlossen werden, dass eine nur 
suboptimale Losung des linearen Optimierungsproblems gefunden wird. Dass es 
sich bei linearen Problemen insbesondere um ein sogenanntes konvexes Problem 
handeln kann, ist prinzipiell bereits bekannt, siehe beispielsweise das Kapitel ..Linear 
Programming" unter dem folgenden Internet-Link der European Business School der 
Universitat Schloss Reichartshausen: 

http://www.ebs. de/LehrstuehleAA/irtschaftsinformatik/NEW/Courses/Semester2/Math 
21. 

Durch die Berucksichtigung der Randwertbedingungen im Rahmen des Begren- 
zungsfaktors wird erreicht, dass als jeweils nachste Iterationslosung wiederum eine 
Losung bestimmt wird, die innerhalb eines zulassigen Wertebereiches liegt, da durch 
den Begrenzungsfaktor die Iterationsschrittweite entsprechend angepasst werden 
kann, so dass eine Randwertbedingung nicht verletzt wird. Durch eine mogliche 
Berucksichtigung der Randwertbedingungen im Rahmen der Bildung des skalierten 
Iterationsgradienten kann insbesondere erreicht werden, dass die Gradientenrich- 
tung moglichst weitgehend so gewahlt wird, dass als jeweils nachste Iterationslosung 
wiederum eine Losung bestimmt wird, die innerhalb eines zulassigen Wertebereiches 
liegt. 

Fur die hier vorliegenden linearen Optimierungsprobleme ist uberdies bekannt, dass 
die optimale, zulassige Losung, die einem Koordinatenpunkt in einem durch 
Randbedingungen begrenzten mehrdimensionalen Raum aller zulassigen Losungen 
entspricht, auf dem Rand dieses begrenzten Raumes liegt. Damit werden der 
skalierte Iterationsgradient und der Begrenzungsfaktor bevorzugt so angepasst, dass 
eine iterative Annaherung an einen optimalen Punkt auf der Begrenzung des 
mehrdimensionalen Raumes der zulassigen Losungen erfolgt. 

Es kann nun insbesondere vorgesehen werden, dass als eine Randwertbedingung 
eine obere Grenze (upper bound) fur eine zulassige Losung definiert wird. 
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Weiterhin kann bevorzugt vorgesehen werden, dass die Bestimmung des 
Iterationsgradienten mit Hilfe einer GaufJ-Elimination erfolgt, was ein sehr schnelles 
Verfahren darstellt. 

Es kann auch insbesondere vorgesehen werden, dass in jedem Iterationsschritt eine 
Skalierung des Iterationsgradienten derart erfolgt, dass eine Gradientenkomponente 
umso kleiner ausfallt, je naher die entsprechende Komponente des Ergebnisses des 
vorherigen Iterationsschrittes einer Randwertbedingung kommt. Es wird damit in 
jedem Iterationsschritt eine neue Skalierung des Iterationsgradienten durchgefuhrt, 
wobei gewisse Komponenten des Gradienten verschwinden, wenn die entsprechen- 
den Komponenten der vorherigen Iterationslosung einer Randwertbedingung sehr 
nahe kommen, d.h. beispielsweise kleiner als ein erstes vordefiniertes AbstandsmafJ 
werden. Dieses erste Abstandsmaft kann auch infinitesimal klein gewahlt werden. 

Weiterhin kann vorgesehen werden, dass die Iterationsphase entweder dann 
beendet wird, sobald das Ergebnis eines Iterationsschrittes mindestens eine Rand- 
wertbedingung uberschreitet und dass das Ergebnis des vorherigen Iterations- 
schrittes als optimaie Losung des Gutekriteriums bestimmt wird. Es wird also die 
Iteration dann abgebrochen, wenn der Algorithmus das Gebiet der zulassigen 
Losungen verlassen wurde, und die letzte zulassige Losung wird als optimaie Losung 
bestimmt. Damit wird auf einfache Weise garantiert, dass in jedem Fall eine 
moglichst optimaie und zugleich zulassige Losung als Endergebnis des Verfahrens 
bestimmt wird. 

Es kann die Iterationsphase aber auch dann beendet werden, sobald das Iterations- 
verfahren gegen eine zulassige Losung konvergiert und sich das Ergebnis eines 
bestimmten Iterationsschrittes von dem Ergebnis eines vorhergehenden Iterations- 
schrittes um weniger als ein zweites vordefiniertes Abstandsmaft unterscheidet, 
wobei das Ergebnis des letzten Iterationsschrittes als optimaie Losung des 
Gutekriteriums bestimmt wird. 
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Nachfolgend wird ein spezielles Ausfuhrungsbeispiel der vorliegenden Erfindung 
erlautert. 



Es wird ein Verfahren zur treibstoffoptimalen Ansteuerung einer Anordnung von 
Dusen an einem Raumfahrzeug betrachtet, welches zur Bahn- und Lageregelung 
des Raumfahrzeuges verwendet wird. Zur Erzeugung von Kraften und Momenten auf 
einem Raumfahrzeug, beispielsweise um simultan Translation und Rotation wahrend 
einer Andockphase oder einer sonstigen Bahn- und Lageregelung regeln zu konnen, 
sind n > 7 Dusen notwendig. Die zugehorigen Ansteuersignale mussen den 
Restriktionen genugen 



positiv zu sein 

kleiner einem Maximalwert (i.a. gleich 1) zu sein 



Daruber hinaus kann bei mehr als 7 Dusen zusatzlich ein Gutekriterium, das i.a. dem 
Treibstoffverbrauch entspricht, optimiert werden. 



Die mathematische Formulierung ftihrt damit auf das folgende lineare Optimierungs- 
problem (LOP): 



Gesucht: Dusenansteuerung a, fur die gilt 



(a) 
(b) 



0 < a <f 
T c a = r 



dim T c = 6 x n, 
(1 ) so dass das folgende Gutekriterium erfullt ist: 
c) Zgj a, = g T a -> min 
wobei 

a: Ansteuervektor 

f: Vektor, der die Maximalwerite der a* enthalt 
g: Vektor der Gewichtsfaktoren 

( f\ f \ 

Tc = jj^jji . Kraft- und Momentenmatrix 

^t c \ y ...,t cn j 

r: Vektor, der die geforderten Krafte und Momente enthalt 
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Zur Anwendung aller LOP Losungsverfahren muss in einer Initialisierungsphase 
zunachst eine zulassige Losung gefunden, d.h. ein Vektor a 2 , der (1a) und (1b) 
erfullt. Mit der sogenannten Singularwertzerlegung (singular value decomposition, 
kurz SVD) von T c 

T c = V 2 Ui T 
mit 

V: ~| 

>- orthogonale Matrizen 
U = (U 1 ,U 2 ): J 

-> T c U 2 = 0 

I = diag <7j > 0: Diagonalmatrix der Singularwerte 

lassen sich alle Losungen von (1b) angeben als 
a 0 = Ui r 1 V T r + U 2 c 

(2) =: U! s + v 

Aus (2) sieht man folgendes: 

(i) zur Realisierung von beliebigen Kraften und Momenten mussen alle o\ grolier 
null sein, d.h. T c vollen Rang haben 

(ii) der erste Summand ist durch r vollstandig bestimmt und stellt die Losung mit 
minimaler Norm von (1b) dar 

(iii) der zweite Summand mit dem noch zu bestimmenden Vektor c dient zur 
Erfullung der Randbedingung (1a) und zur Minimierung von (1c) 
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Aus der Tatsache, dass mit dem Diisensatz sowohl positive wie negative r 
realisierbar sein mussen, folgt aus (2), dass es c-\ und c 2 geben muss, so dass 

(a) ai := Ui si + U 2 Ci > 0 

(b) a 2 := Ui (-Si) + U 2 c 2 > 0 

(3) 

(c) -> U 2 (C! + c 2 ) =: U 2 c p > 0 

d. h. die Existenz beliebig vieler Vektoren c p mit 

(4) v p = U 2 c p > 0 

gesichert ist. Nach Wahl eines bestimmten (a priori fixierbaren) v p la&t sich 
somit ao positiv machen gemali 



(5) 



(a) a 0 = Ui s + v p ki > 0, 

(b) k n = max LHi*)l±± f g „ q q 4 > 0 

vp, 



wobei e aus numerischen Grunden fur die Anwendung der nachfolgenden 
Optimierungsschritte eingefOhrt wurde. 

Fur grofie rechte Seiten r kann es vorkommen, dass (1a, b) keine Losung a < f hat, 
es wird daher das um x s erweiterte Problem betrachtet 



(a) T c a = r (1 - x s ), 



(b) 



0 < a <f 
0 < x s < 1 



(6) (c) ^ (Tc, r) 



= r 



(d) g T a + g s x s -> min 

(e) gs > g T f, damit x s zu null wird, falls r nicht zu graft ist 
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Damit lasst sich nun auch die obere Begrenzung einhalten und der gesuchte 
zulassige Startwert fur a berechnen zu 



(a) a z = a 0 (1 - x sz ) = Ui s (1 - x sz ) + v p ki (1 - x S2 ) 
f w, w > 0 
0, w < 0 



(7) (b) x sz = 



( C ) w = max \ 1 - — — - 1 

e stellt hier ein erstes, beispielsweise infinitesimales Abstandsmafc dar. 

Alle folgenden Betrachtungen beziehen sich auf das erweiterte System (6), wobei 
aber die ursprungliche Bezeichnung nach (1) der Einfachheit halber beibehalten wird. 

Zur Losung des LOPs erfolgt nun ein zweiter Verfahrensschritt, namlich eine 
Optimierung des Gutekriteriums (1c) bzw. (6d), die iterativ wie folgt durchgefuhrt 
wird: 



a i+1 = aj - D, U 2 (i) U 2 (I)T D, g ■ k 
(a) ' 

- B\ - Vgj 



y 

Vni ■ k 



Dabei stellt v gi den Iterationsgradienten dar, der in jedem Iterationsschritt 
skaliert wird, d.h. es wird in jedem Iterationsschritt die Gradientenrichtung neu 
bestimmt. Weiterhin stellt k einen Begrenzungsfaktor fur die Iterationsschritt- 
weite dar, der wie folgt bestimmt wird: 



(8) 



(b) k = min (k Uf k 0 ) (1 - e) 



(c) 



k n - min 

vQ)<o 



U) 



fj 



& J 



k., = min 

,W>0 



^ I 

V & J 
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Diese Wahl von k unter Berucksichtigung der Randwertbedingung 0 < a < f 
sichert, dass a i+ i zulassig bleibt. 













[ fj) 





(e) U 2 (i) : Nullraumvektoren von T c Dj = V (i) U^ )T 

ai = a z , berechnet nach (7) 

Die wesentliche Idee in (8) ist die standige Skalierung des Problems mit D\ und das 
anschlieftende Fortschreiten in die so modifizierte, auf U 2 {,) projizierte negative 
Gradientenrichtung, womit aufgrund der bekannten Struktur des Problems als 
konvexes lineares Optimierungsproblem mit einer optimalen Losung auf dem Rand 
garantiert wird, dass das Gutekriterium in jedem Iterationsschritt reduziert wird. Die 
Iteration wird bevorzugt dann abgebrochen, wenn der Betrag von v gi k eine 
vorgegebene Schranke als zweites Abstandsmafi unterschreitet, d.h. aj sich kaum 
noch andert. 



Eine besondere Erweiterung des vorliegenden Verfahrens gegenuber dem Vorgehen 
bei Karmakar besteht darin, dass ein zusatzliches Randwertproblem mit in jedem 
Iterationsschritt berucksichtigt wird, hier das Einbeziehen der oberen Grenze f (upper 
bound problem) durch Hinzunahme des 2. Faktors in Dj und Berucksichtigung der 
oberen Grenze f im Rahmen des Termes k 0 in der Berechnung von k nach (8b). 
Bislang wurden bei Karmakar-Verfahren ublicherweise komplexe Erweiterungen des 
linearen Optimierungsproblems um Schlupfvariablen vorgesehen mit dem Nachteil, 
dass sich dann die Dimension des zu losenden Problems deutlich erhoht. Hier stellt 
das vorliegende Verfahren eine wesentliche Vereinfachung dar. Aufierdem entfallt 
die bei Karmakar sehr aufwandige Bestimmung einer zulassigen Losung in der 
Initialisierungsphase durch die hier vorgeschlagene, dem hier vorliegenden Problem 
besser angepasste Initialisierungsphase. 
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Ein weiterer vorteilhafter Verfahrensschritt des hier beschriebenen Verfahrens liegt 
also in der laufenden Berechnung der v gi , die bevorzugt nicht uber die SVD von T c Di 
erfolgt, sondern wegen (benutzt wird hier die folgende, vereinfachte Darstellung: T c 
fur T c Di, D fOr Di) 

(a) M = T C T C T = V E 2 V T 

(b) U 1 U 1 T = T C T M- 1 T, 



v g = D U 2 U 2 T Dg = D(l - UiUiJ^Dp 
(9) = D ( I - T C T M" 1 T c ) b b 



(c) 
mit 



= D (b - T C T x) 



(d) Mx = T c b 

uber die Losung von (9d) fur x mittels einer GauR-EHmination erfolgt. Dieses 
Verfahren ist deutlich schneller als eine SVD. Ebenso wird fur a 0 auch LM nicht 
ermittelt, sondern U-iS direkt berechnet gemaB 

(e) UiS = Ui E" 1 V T r = T C T M" 1 r 

Somit besteht das Verfahren in der Ausfuhrung der Rechenschritte 



Gig (5) mit (9e) 
Gig (7) 



Initialisierung 
Gig (8) mit (9c, d) } Iteration 



Abschliefcend sind weitere Vorteile des vorgeschlagenen Verfahrens gegenOber 
einem Simplex-Verfahren nach dem Stand der Technik zusammengestellt: 

eine Berechnung einer zulassigen Losung kann mit minimalem Aufwand 
erfolgen, (sogar a z kdnnte bereits als Dusenansteuerung benutzt werden, 
ware allerdings nicht treibstoffoptimal). 
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die anschliefcende Optimierung kann gegebenenfalls auf wenige Schritte 
begrenzt werden, wenn z.B. die Leistung des Bordrechners begrenzt ist. 

die Erfahrungen zeigen, dass das Optimum im Mittel mit erheblich weniger 
Rechenoperationen erreicht wird und generell weniger Speicheraufwand 
benotigt als bei einem Simplex-Verfahren. 

bei einem Dusenausfall ist Tc lediglich urn die entsprechende Spalte zu 
reduzieren und ein entsprechender Vektor v P abzuspeichern, d.h. es ist in 
diesem Fall nur ein minimaler zusatzlicher Speicherbedarf erforderlich, im 
Gegensatz zum Stand der Technik. 

durch die Einbeziehung des upper bound problems werden in der Praxis 
Probleme besonders bei begrenzter Dusenstellkapazitat vermieden. 
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Patentanspriiche 

1. Verfahren zur treibstoffoptimaien Ansteuerung einer Anordnung von Dusen an 
einem Raumfahrzeug, wobei zum Auffinden einer treibstoffoptimaien Losung fur 
die Ansteuerung die Losung eines linearen Optimierungsproblems erfolgt mit 

einer Initialisierungsphase zum Auffinden einer ersten zulassigen Losung des 
linearen Optimierungsproblems und 

einer anschliefienden Iterationsphase, in der ausgehend von der zulassigen 
Losung des linearen Optimierungsproblems eine iterative Optimierung eines 
Gutekriteriums erfolgt, 

dadurch gekennzeichnet, dass 

bei jedem Iterationsschritt die Bildung eines skalierten Iterationsgradienten 
erfolgt und 

der Iterationsgradient mit einem Begrenzungsfaktor fur eine maximale Iterati- 
onsschrittweite multipiiziert wird, welcher unter Berucksichtigung mindestens 
einer Randwertbedingung fur eine zulassige Losung gebildet wird. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass als eine Randwert- 
bedingung eine obere Grenze fur eine zulassige Losung definiert wird. 

3. Verfahren nach einem der Anspruche 1 oder 2, dadurch gekennzeichnet, dass 
die Bestimmung des Iterationsgradienten mit Hilfe einer Gaufi-Elimination 
erfolgt. 

4. Verfahren nach einem der Anspruche 1 bis 3, dadurch gekennzeichnet, dass in 
jedem Iterationsschritt eine Skalierung des Iterationsgradienten derart erfolgt, 
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dass eine Gradientenkomponente umso kleiner ausfallt, je naher die entspre- 
chende Komponente des Ergebnisses des vorherigen Iterationsschrittes einer 
Randwertbedingung kommt. 

5. Verfahren nach einem der Anspruche 1 bis 4, dadurch gekennzeichnet, dass 
die Iterationsphase beendet wird, sobald das Ergebnis eines Iterationsschrittes 
mindestens eine Randwertbedingung uberschreitet und dass das Ergebnis des 
vorherigen Iterationsschrittes als optimale Losung des Gutekriteriums bestimmt 
wird. 

6. Verfahren nach einem der Anspruche 1 bis 4, dadurch gekennzeichnet, dass 
die Iterationsphase beendet wird, sobald das Iterationsverfahren gegen eine 
zulassige Losung konvergiert und sich das Ergebnis eines bestimmten 
Iterationsschrittes von dem Ergebnis eines vorhergehenden Iterationsschrittes 
um weniger als ein vordefiniertes Abstandsmafi unterscheidet, wobei das 
Ergebnis des letzten Iterationsschrittes als optimale Losung des Gutekriteriums 
bestimmt wird. 
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Zusammenfassung 

Verfahren zur treibstoffoptimalen Ansteuerung einer Dusenanordnung 

Beschrieben wird ein Verfahren zur treibstoffoptimalen Ansteuerung einer Anordnung 
von Dusen an einem Raumfahrzeug unter Losung eines linearen Optimierungs- 
problems mit einer Initialisierungsphase zum Auffinden einer ersten zulassigen 
Losung und einer anschliefienden Iterationsphase, in der ausgehend von der 
zulassigen Losung eine iterative Optimierung eines Gutekriteriums erfolgt. Bei jedem 
Iterationsschritt erfolgt die Bildung eines skalierten Iterationsgradienten und der 
Iterationsgradient wird mit einem Begrenzungsfaktor fur eine maximale Iterations- 
schrittweite multipliziert, welcher unter Berucksichtigung mindestens einer Randwert- 
bedingung fur eine zulassige Losung gebildet wird. 
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